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ABSTRACT 

We introduce a novel theoretical model to explain the long-standing puzzle of the 
nature of the microlensing events reported towards the Large Magellanic Cloud (LMC) 
by the MACHO and OGLE collaborations. We propose that a population of tidally 
stripped stars from the Small Magellanic Cloud (SMC) located ~4-10 kpc behind a 
lensing population of LMC disk stars can naturally explain the observed event dura- 
tions (17-71 days), event frequencies and spatial distribution of the reported events. 
Differences in the event frequencies reported by the OGLE (~0.33 yr _1 ) and MACHO 
(^1.75 yr _1 ) surveys appear to be naturally accounted for by their different detection 
efficiencies and sensitivity to faint sources. The presented models of the Magellanic 
System were constructed without prior consideration of the microlensing implications. 
These results favor a scenario for the interaction history of the Magellanic Clouds, 
wherein the Clouds are on their first infall towards the Milky Way and the SMC has 
recently collided with the LMC 100-300 Myr ago, leading to a large number of faint 
sources distributed non-uniformly behind the LMC disk. In contrast to self-lensing 
models, microlensing events are also expected to occur in fields off the LMC's stellar 
bar since the stellar debris is not expected to be concentrated in the bar region. This 
scenario leads to a number of observational tests: the sources are low-metallicity SMC 
stars, they exhibit high velocities relative to LMC disk stars that may be detectable 
via proper motion studies, and, most notably, there should exist a stellar counterpart 
to the gaseous Magellanic Stream and Bridge with a V-band surface brightness of 
> 34 mag/arcsec 2 . In particular, the stellar Bridge should contain enough RR Lyrae 
stars to be detected by the ongoing OGLE survey of this region. 

Key words: gravitational lensing — dark matter — Galaxy: halo — Galaxy: struc- 
ture — Magellanic Clouds — galaxies: interactions 



1 INTRODUCTION 

IPaczvnskH (|l986l ) suggested that massive compact halo ob- 
jects (MACHOs) could be found by monitoring the bright- 
ness of several million stars in the Magellanic Clouds (MCs) 
in search for microlensing by unseen foreground lenses. De- 
spite the seemingly daunting scale of such a project, there 
have been a number of surveys conducted towards the MCs 
in se arch of microlensin g events, such as the MACHO sur- 
vey (|Alcock et al.|[2000l ). the Experi ence pour la Recherche 
D'Objets Sombres survey (EROS; iTisserand et~aH 120071 1 
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and the Optical Gr avitational Lensing Experiment (OGLE; 
lUdalski et al1ll997n . The goal of these surveys was to test 
the hypothesis that MACHOs could be a major component 
of the dark matter halo of the Milky Way (MW). After 20 
years of work, a number of microlensing events have been 
detected towards the Large Magellanic Cloud (LMC), but a 
clear explanation for their origin remains elusive. 

The results from the MACHO, OGLE and EROS sur- 
veys (summarized in § 12. 2[) have demonstrated that objects 
of astrophysical mass (10 -7 - 10 M B ) do not comprise the 
bulk of the dark matter in the halo (|Alcock|[200alMoniezl 
1201(f) . 

The idea that the MW's dark matter halo is completely 
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populated by MACHOs has been ruled out based on el- 
emental abundances and constraints on the total baryon 
fraction from baryon acoustic oscillations observed in the 
cosmic microwave background (|Komatsu et al.ll201ll ). If the 
observed events result from lensing by MACHOs, the im- 
plied MACHO mass fraction (assuming these objects are 
^0.5 Mq in m ass) of the MW dark matter halo is ~16 % 
jBennettll2005h . The nature of such MACHOs is unknown. 
Low mass hydrogen burning stars have been ruled out as a 
possib ility by red star c ounts from the Hubble Space Tele- 
scope (|Flvnn et alJll996h . White dwarfs have also been pro- 
posed as potential candidates, supported by claims that > 
2 % of the halo is possib ly made up of old white dwarfs 
|Oppenheimer et alj 2001). However, it appears likely that 
the obse rved white dwarfs are part of the thick disk of 
the M W jTorres et al.ll2002l;lFlvnn et al.ll2003l ; ISpagna et all 



|2004 etc.). iBrook et ahl (|2003l ) also rule out the possibility 



of a white dwarf halo fraction in excess of 5% based on the 
inability of an initial mass function biased heavily toward 
white dwarf precursors to reproduce the observed C/N/O 
abun dances in old MW stars (see also iGibson fc Mould! 
1 19971 ). 

Alternatively, the observed microlensing signal might 
be a result of lensing by normal stellar populations. Self- 



detected by the OGLE team (LCalchi Novati et al. 


I200S ; 


Calchi Novati & Mancinil l201ll; IWvrzvkowski et al. 


2011b: 


Sahu 1994). Also binary star lenses were found to be associ- 


ated with the LMC/SMC (iBennett et al.lll99d:lAfonso et al.1 


200d: Udalski et al. 19941; Dominik & Hirshfeldl 


1996; 


Mao & Di Stefand 199a Alard et al. 19951). SahrJ 


(00031) 



raise a number of additional arguments disfavoring a halo 
origin for the lensing events. In particular, if the lenses 
are MACHOs in the halo, the durations of the events in 
the LMC should be similar to those in the SMC. How- 
ever, the SMC event durations (75-125 days) are observed 
to be much longer than those towards the LMC (17-71 
days). This discrepancy can be reconciled in models where 
the lenses originate within the Clouds themselves since the 
SMC has a much larger line-of-sight depth than the LMC 
(lHaschke et al.l [20121 : ISubramanian fc Subramaniaml 120121 : 
ICrowl et al.ll200ll ). 

However, attempts to explain the MACHO team's 
larger number of reported events using a variety of disk and 
spheroid models for the LMC and MW have failed to explain 
their reported mic r olensing optica l depth at the 99.9% level 
llGvuk et al.ll200d ; IBennettl [20051 . and references therein). 



I Gould ( 1995 ) presents simple dynamical arguments against 
sel f-lensing models by L MC stars (such as those proposed 
bv lSahu|[l99l : IWulll994h . finding that the required velocity 
dispersions would be much larger than observed for known 
LMC populations in order to reproduce the observed lensing 
optical depth. 

Non-standard models for the LMC's structure have 
been explored to tr y and get around such constraints. 
I Zhao fc Evansl |200d ) suggest that the off-center bar in the 
LMC is an unvirialized structure, misaligned with and off- 
set from the plane of the LMC disk by 25 degrees. They 
argue that this offset allows the bar structure to lens 
stars in the LMC's disk. As discussed in iBesla et al.l ^Old - 
hereafter, B12 ) , our simulations agree with the theory of 
I Zhao fc Evansl (2000) that the bar of the LMC has become 



warped by its recent interactions with the SMC - but the 
simulated warp of ~10 degrees is less severe than what they 
require. Moreover, microlensing events have been identified 
off the bar of the LMC, meaning that an offset bar cannot 

be t he main explan a tion. 

iMancini et aU [[20041 ) refine the iGouldl (|l995l ) ar- 
gument and account for our updat ed understand i ng of 

" [200l|; 
they 



the geometry of the L MC's disk l|van der M arc 



ir M arel 

Ivan der Marel et alJlioolj ). Like lZhao fc Evansl l|200fj| ) 

argue that a non-coplanar bar strengthens the self-lensing 
argument. But ultimately they find that self-lensing cannot 
account for all the events repo rted by the MACHO team 
|Calchi Novati fc Mancinill201ll). 



Weinberg! l|200d ) and lEvans fc Ker ins (2000) claim that 



the LMC may have a non-virialized stellar halo, differing 
from a traditional stellar halo in that it has a low velocity 
dispersion (for that reason it is referred to as a "shroud" 
rather than a halo). A stellar halo is reported to have 
been detected around the LMC, extendin g as much as 20 
degrees (20 kpc) from the LMC center (Maicwsk i et al.1 
20091: iMunoz et al.l[200d . and Nidever et al. 2012, in prep.fl 
Peicha fc StanekM 2009) used RR Lyrae stars from the 



OGLE III catalogue to claim that such an LMC stellar halo 
is structured as a triaxial ellipsoid elongated along the line of 
sight to the observer, which would aid in explaining the mi- 
crolensing events. But the total mass of such a halo is likely 
no m ore than 5% of the mass of the LMC (|Kinman et alj 
Il99ll) . and is thus unlikely to explain the o bserved microlens- 
ing optical depth: lAlves fc Nelson! (|200d ) place a limit on 
the fractional contribution to the optical depth from such a 
low m ass LMC stellar halo component to 20%. [Cvuk et all 
(2000) argue that in order to explain the total optical depth, 
the mass of the stellar halo must be comparable to that of 
the LMC disk + bar, a scenario t hat is ruled out by obser- 
vati ons (see also the discu ssion in lAlcock et al.l feoOl), 

IZaritskv fc Lml (|l997f ) claim to have detected a fore- 
ground intervening stellar population with a stellar surface 
density of E ~ 16 Mq pc~ 2 , w hich could serve as the lensing 
stars. However, Gould ( 19981 ) argues that such a population 
would need an anomalously high mass-to-light ratio in or- 
der to explain a large fraction of the observed events and 
remain unseen. IBennettl ([ 1998h also point out that a fore- 
ground stellar population that could explain the LMC opti- 
cal depHijwojiklreqjrire more red clump star s than inferred 
bv lZaritskv fc Lin! j| 19971 ) . Furt hermore, both lGallartl (| 19981) 



and 



Beaulieu fc Sackett (| 19981 ) argue that the feature seen 



by IZaritskv fc Lid (|l997T l in the color magnitude diagram 



is a natural feature of the red giant branch, rather than an 
indication of a foreground population. 

Thus, 10 years after the MACHO team's findings, we 
are still left with unsatisfactory solutions to the nature of 
the observed microlensing events. 

In this study, we consider the possibility that LMC 
stars are not the sources of the microlensing events, but 
rather the lenses of background stellar debris stripped from 
the SMC. This theoretical model is based on our recently 



1 Although it is un clear whether the ve ry extended stellar pop- 
ulation detected bv lMunoz et al. (2006) includes main sequence 
stars (Olszewski et al. 2012 in prep), shedding doubt on their 
membership as part of an extended LMC stellar halo. 
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published simulations of the interaction history of the MCs 
and the formation of the Magellanic Stream, a stream of HI 
gas that trails behind the MCs 150 degrees across the sky 
jNidever et al.ll2010l ; iPutman et"al]|2003l : iMathewson et all 
1 1974 ), in a scenar io where the MCs are o n their first infall 
tow ards o u r MW l|Besla et al.ll20ll I2012T ). 

IZhad l| 19981 ) (hereafter Z98) suggested that the mi- 
crolensing events towards the LMC could be explained by 
tidal debris enshrouding the MCs. This debris is assumed to 
be associated with the formation of the Magellanic Stream 
and could be either in front or behind the LMC, i.e. acting 
as sources or lenses. In particular, Z98 was able to predict a 
lower concentration of events towards the LMC center than 
the self-lensing models, which rely strongly on the bar. Our 
model differs in that the debris (source) stars are a combina- 
tion of SMC stars captured by the LMC over previous close 
encounters between the MCs (<1 Gyr ago) and a population 
of younger stars pulled out from the SMC (or forming in the 
Bridge) during the rec ent collision with the LMC 100-300 
Myr ago (Model 2 from lBesla et alj|2012ft . Z98, on the other 
hand, hypothesize that the debris is a significantly older stel- 
lar population associated with the Magellanic Stream, span- 
ning a great circle about the MW. The resulting kinematics 
and spatial distribution of the sources in our model is con- 
sequently quite different than in Z98. Moreover, the great 
circle picture of stellar debris is inconsistent with a first in- 
fall scenario, as the MCs have not yet completed an orbit 
about the MW. 

In the B12 model, the Magellanic Stream owes its 
origin to the action of LMC tides stripping material 
from the SMC before the MCs have been accreted by 
the MW. This scenario is in sharp contrast to prevail- 
ing models for the origin of the Stream, which all rely 
to some degree on the action of MW tides at previous 
pericentric passages to explain its full extent and the ex- 
istence of material leadin g to the MCs, referred to as 
the L e ading Arm Featur e (Heller & Rohlfs 1994 ; [Lin et alj 



19951; IGardiner fc Noguchil Il996l; iBekki fc Chibal 120051; 
Connors et al.ll2 005: Mastropi etro et al.l 2005; Ruzicka et alj 
20ld ; iDiaz fc Bekkil l201ll . I2012T I. Such a scenario is at 



odds with a first infall, which is the expected orbital so- 
lution for the MCs based on updated c osmological mod- 
els for the dark matter ha l o of the MW dBesla et all 120071; 
iBovlan-Kolchin et al.ll201ll ; lBusha et al.l201ll) and new HST 
proper motion measurements of the MCs ( Kallivavalil et all 
l2006al lbl. Kallivayalil et al., 2012 in prep). 

The B12 LMC-SMC tidal model should remove stars as 
well as gas; however, to date no old stars have been detected 
in the Stream or the Magellanic Bridge that connects the two 
galaxies. Indeed the apparent absence of stars has been the 
main argument in favor of hydrodynamic models, where the 
Stream forms from ram pressure stripping by the MWs dif- 
fuse hot gaseous halo or from stellar outflows. All tidal mod- 
els for the Magella nic Stream predict a stellar counterpart to 
some degree fe.g., |Piaz fc Bekkill2012l ; IGardiner fc Noguchil 
1 19961 ). whereas the hydrodynamic models naturally explain 
their absence (e.g., Nidever et al.l 120081 ; iMastropietro et alj 
2005; M oore fc Davisl [19941 ) . In this work we demonstrate 
that the non-detection of a stellar counterpart to the Mag- 
ellanic Stream is reconcilable with the B12 model, as tides 
remove material from the outskirts of the original SMC's 
extended gaseous disk, where the stellar density is low. 



The corresponding surface brightness of the predicted stellar 
counterpart to the Magellanic Stream and Bridge is below 
the sensitivity of existing surveys, but should be observ- 
able by future and ongoing photometric surveys. In fact, 
such a debris field may already have been detected: re- 
cently, lOlsen et alj (|201ll ) discovered a population of metal 
poor RGB stars in the LMC field that have different kine- 
matics from tho se of local stars in the LMC disk (see also 
iGraff et al.ll2000h . In B12 we showed that such kinematically 
distinct debris is expected from tidal LMC-SMC interactions 
and represents ~1.5% of the LMC's disk mass. 

The overall goal of the work presented here is to use 
the simulated stellar debris predicted to exist behind the 
LMC by the B12 models to determine the corresponding 
microlensing event properties (durations, frequency, distri- 
bution) towards the LMC owing to the lensing of this debris 
by LMC disk stars. These values will be compared directly 
to the MACHO and OGLE survey results. 

We begin our analysis with a description of the results of 
the three main microlensing surveys (§ [2} and our method- 
ology and simulations (§|3J. We follow with a detailed char- 
acterization of the surface densities and distribution of the 
sources and lenses predicted by the B12 tidal models for the 
past interaction history of the MCs (§|3J- We will henceforth 
refer to SMC debris stars as sources and LMC disk stars as 
lenses. We then estimate the expected effective distance be- 
tween the simulated lens and sources (§[STTJ and the relative 
velocities transverse to our line of sight f§ 15. 2[) . Using these 
quantities we derive the event durations (§ 15 . 3f) and event 
frequency (§ I5.4f) from lensed SMC debris and compare di- 
rectly to the MACHO and OGLE results described in § l2~2l 
Finally we outline a number of observationally testable con- 
sequences of the presented scenario (§ ED • 



2 MICROLENSING OBSERVATIONS 

First we provide an overview of the relevant equations (§ I2.1[1 
and summarize the results of the MACHO, OGLE and 
EROS surveys f§ 12. 2p . In this study we aim to compute the 
duration and number of expected microlensing events that 
would be observable over the duration of the MACHO and 
OGLE surveys from two different models of the interaction 
history of the LMC-SMC binary system, as introduced in 
the B12 study. 



2.1 Microlensing Equations 

The term microlensing refers to the action of compact low 
mass objects along the line of sight to distant sources, where 
"compact" refers to objects smaller than the size of their own 
Einstein Radii, and "low mass" refers to objects between 1 
M® and 10 3 Mq . The resulting splitting of the image of the 
distant source is typically unresolvable, appearing instead as 
an overall increase in the apparent brightness of the source. 

The microlensing optical depth, r, is defined to be the 
instantaneous probability that a random star is magnified 
by a lens by more than a factor of 1.34 (corresponding to an 
impact parameter equal to the Einstein radius). The opti- 
cal depth depends only on the density profile of lenses and 
the effective sep aration between the lens and sources (D) 
l|Paczvnskilll986l ): 
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D 



DlDsl 
D s 



(1) 



where Dl is the distance between the observer and the lens, 
Ds is the distance between the observer and the source and 
Dsl = Ds — Dl- Throughout this study we compute quan- 
tities in gridded regions spanning the face of the LMC disk, 
as defined in § 14.21 We take Dl as the mean distance to 
LMC disk stars in each grid cell but compute Ds explicitly 
for each SMC star particle individually. We identify source 
debris stars as those SMC stars with Dsl > to ensure that 
we only choose source stars that are behind the LMC disk 
on average. 

The microlensing optical depth is formally defined as: 



4ttG 



t = 6.0 x 10" 



D 
kpc 



M /pc 2 



(2) 



(3) 



The optical depth reflects the fraction of source stars that 
are magnified by a factor A > 1.3416 at any given time. 
Observationally, r is estimated as: 



fobs 



1 7T 

E2 



E 



tj 



(4) 



where E = 6.12 x 10 7 object years is the total exposure 
time (10.7 x 10 6 objects times the survey duration of 2087 
days), ti is the Einstein ring diameter crossing time of the 
ith event, and £(fj) is its detection efficiency, which peaks at 
0.5 for durations of 200 days for the MACHO survey. The 
value of E is dependent on the assumed source population. 
Thus, r can be approximated as: 



2A S 



(5) 



source 



where T is the event frequency and t e is the event duration. 
We aim to derive these two quantities from the simulations. 

From equation ([5]), it is clear that the value of r is 
strongly dependent on the adopted definition of the source 
population. In the model explored here, we assume the 
source population consists of stars stripped from the SMC. 
This differs markedly from the choices made by the OGLE, 
MACHO and EROS teams, who assumed that the source 
population consists of LMC stars. As such, our value of 
-/Vsource is significantly smaller than that adopted by these 
surveys. Consequently, our computed lensing probability, r, 
cannot be directly compared to the expectations from these 
surveys. Instead we must compare our simulations to the 
quantities that were observed directly, i.e. the event dura- 
tions and event frequency. 

The event duration (t e ) is defined as the Einstein radius 
crossing time: 



U = 



Re 



(6) 



After substituting for the Einstein radius R e , we get 
(D_\° B ( V 



t e = 4944.7 



\kpc J ykm/s 



Ml ens 



M 







days. (7) 



Following lAlcock et all l|200Ch . we do not include event 
durations longer than 300 days. Note that the MACHO 
team adopt a non-standard definition for the event dura- 
tion, which is twice the value we adopt here. We compute 
t e using the average value of (M^ s ) from an adopted initial 
mass function, as described later in the text. 

To compute t e we need to substitute both D (as defined 
in equation [1]) and V±, the relative velocity of the source and 
lens perpendicular to our li ne-of-sight and proje cted in the 
plane of the lens. Following lllan fc Gouldl (|l996h . 



V± = (V L - Vo) - (V s 



Vo) dS> 



(8) 



where Vs, Vl and Vo are the velocities of the source, lens 
and observer, respectively, perpendicular to our line-of-sight. 
These are defined by taking the cross product of the Galacto- 
centric velocity of each particle with the normalized line-of- 
sight position vector of that particle. Vo is defined by taking 
the cross product of the observer velocity, Vo = Vlsr + V pec , 
with the average normalized line of sight position vector of 
the LMC particles. 

Recently, the values of both Vlsr and Vpec have come 
unde r debate. The standard IAU value l|Kerr fc Lvnden-Belll 
Il986h for the circular velocity of the Local Standard of 
Rest (LSR) is Vlsr = 220 km/s . However, models based 
on the proper motion of Sgr A* (|Reid fc Brunthalerll2004h 
and m asers in high-mass star-formation regions (|Reid et alj 
2009) have suggested t hat the circular velocity may be 
higher. iMcMillarj (|201lh has presented an analysis that in- 
cludes all relevant observational constraints, from which he 
derived Vlsr = 239±5 km/s . In the following, we adopt this 
value for Vlsr and his distance for the Sun from the MW 
center of R© = 8.29 ± 0.16 kpc. For the peculiar velocity of 
the S un with respect t o the LSR we take the recent estimate 
from lSchonrich et all l|201Ch : V pec = (U pec , V pcc , W pC c ) = 
(11.1, 12.24, 7.25), with uncertainties of (1.23, 2.05, 0.62) 
km/s. 

This was not the approach adopted in B12. The revised 
Vlsr and Vpec imply a solar velocity in the Galactocentric- 
Y direction that is 26 km/s larger than that used b y 
iKallivavalil et al.l (|2006aT ) and IKallivavalil et all (|2006bh . 
This directly impacts the dete rmination of the 3D v eloc- 
ity vector of the LMC (see also. IShattow fc Loebll2009l ). We 
discuss the consequences of these revised values for the sim- 
ulations analyzed in this study in § [3] 

Because stellar particles representing source/lens stars 
will exhibit a wide variation in V± and D, we also expect 
to find a range of plausible t e in a given field of view. This 
range of values can be compared directly to the MACHO 
and OGLE results. We compute the average value of t e by 
integrating over the probability distribution of P(D 5 /Vj_): 



(t e ) = 4944.7( 



M,: 



/ 



V, 



\kpc / \km/s 



days 



P(D°- 5 /Vx)d(D°- 5 /V±). 



(9) 



/ \0.5 

To compute ( I N 1 ™ B J ) , we define a mass spectrum 
following iBastian et al? l|2010r i and ISumi et ail l|201ll . sup- 
plementary information) for the low mass end: 
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dN 



where 



a 
a 
a 



0.3 (0.01 M Q < M, < 0.075 M ) 
1.3 (0.075 M Q < M, < 0.5 M ) 
2.35 (0.5 M Q < M, < 10 M ) . 



oc M- a (10) 

(11) 
(12) 

(13) 

This IMF follo ws a Salpeter i ndex (a - 1 = 1.35) at 
the high mass end (Salpeter 1955) and a shallower power 
law (a — 1 = 0) for lower mass stars down to the hydrogen 
burning limit of 0.075 Mq. 

We adopt a mass range of 0.01-3 Mq for lensing stars. 
Stars more massive than 3 Mq are bright enough to have 
been easily observable. The lower limit accounts for the pos- 
sibility that brown dwarfs could act as lenses; brown dwarfs 
are as common as main sequence stars in the Milky Way, 
meaning that a large fraction o f Galactic bulge microlens- 
ing events owe to brown dwarfs l|Sumi et al.ll20lj) . Reliable 
statistics for the mass function of b rown dwarfs exist down 
to about 0.01 M (|Sumi et al.ll201lf ). With these given mass 
limits, the average square root of the lensing mass is: 



J dm 

f^dm 

J dm 



0.61 M Q - 5 



(14) 



Finally, the event frequency, F, is computed from equa- 
tions ©, © and (0 as : 



M /pc 5 



M k 



(15) 

where Ei ens is the average surface mass density of lenses in 
each grid. To determine JV SOU rce, the number of source stars 
per grid cell, we normalize the IMF given in eqn [10] using 
the total mass of stellar debris from the SMC per field of 
view. We then compute the number stars within a subset 
mass range of 0.075-3 M , which we adopt for source stars. 
The lower limit corresponds to the hydrogen burning limit 
and the upper limit to ~ the mass of a main sequence A 
star. 

To compute the average number of events expected per 
year within a given field of view, we must account for the 
probability distribution of P(Vj_D - 5 ) for the stripped SMC 
stellar particles (sources) within that region: 



(T) = 2.8 x lO^AT, 



source 



M /pc 



Mr. 



J V±D - 5 P(V±D - 5 )d{V ± D°- 5 )yT~ 1 . 



(16) 



Here, P(Vj_D ' 5 ) is the normalized probability distribution 
for the product of V±D 5 using equations ([TJ and © for 
all sources in a given field. This method assumes that all 
sources are screened by the same lens population, meaning 
that Ei cns is assumed to be uniform across the field of view 
(grid cell) considered. 

Again, using the IMF in eqn 1101 the relevant average 
lensing mass (mass range 0.01-3 M ) in eqn ll6l is: 



(Mens 5 ) 



[ m -0.5dN dm 
J dm 

f^dm 

J dm 



= 2.11M 







(17) 



2.2 Results of the MACHO, EROS and OGLE 

Surveys 

The MACHO survey examined the central ~14 deg 2 region 
of the LMC over 5.7 years and identified 17 c andidate mi- 
crolensing events (set B of lAlcock et aljfeooch . resulting in 
an estimate for the microlensing optical depth towards the 
LMC of t = 1.21q 3 x 10~ 7 . Since the original publication, a 
number of these candidates have been rul ed out as tran - 
sient events, such as SNe (e. g. Event 2 2; Bennettl |2005| ). 



variable stars (e.g. Event 23: [Tisserand! [20051 ') or as lens- 
ing events by st ars in the thick disk of the MW (e.g. Event 
20 and Event 5: iKalliyavalil et al.ll2006l: iNguven et aT]l2004l : 
IPopowski et alil2003l : brake et al.ll2004h . OGLE III data re- 
cently shows that Event 7 exhibited a few additional bright- 
ening episo des, excluding it from bein g a genuine microlens- 
ing event l|Wvrzvkowski et al.l 1201 laT). Som e authors have 
made different cuts: iBelokurov et al.l (|2004l ) claim that only 
7 events are real. 

After the listed corrections, 10 of the MACHO events 
survive as likely candidates and are listed in Table [1] along 
with the duration of each event and source magnitudes. The 
e vent number is in re ference to the event numbers assigned 
m lAlcock et al1 |2000). The corresponding event frequency is 
10/5.7 ~1.75 per year. This is lower than would be expected 
if the dark halo of the MW was mad e primarily of objects 
with masses of ~ 5M (|Alcockll20o"9T ). 

The EROS-1 survey covered a 27 deg 2 region contain- 
ing the LMC bar and one field toward the S MC for 3 years. 
Two events were reported towards the L MC (lAubourg et alj 



1 1993h . but both have been rejected llAnsari et al.l 119951: 
[ Renault et alj 1 19971 : ILasserre et alj l2000l : Tisserand et alj 
l2007ft . The followup EROS-2 survey covered a larger region 
of the LMC, spanning 84 deg 2 , over 6.7 years found no candi- 
date microlensing event towards the LMC (|Tisserand et alj 
120071 ): combined with the EROS-1 results, this implied a 
microlensing optical depth of less t han r < 0.36 x 10 -7 
|Tisserand et al.l 120071 : iMoniedlioioh . The discrepancy be- 
tween the EROS-2 and MACHO results likely owes to the 
EROS team's choice to limit their candidate sources to only 
clearly identified bright stars in sparse fields throughout the 
LMC (0.7 x 10 7 stars over 84 deg 2 ), rather than including 
blended sources and faints s tars in dense fields (MACHO: 
1.2 x 1 7 stars over 1 4 deg 2 ) (jMonicz 2010). As also pointed 
out bv lMoniezI d2010t). on l y 2 of the 17 MACHO candidates 
in set B of Alcock et alj l|2000l ) were sufficiently bright to 



be compared to the EROS sample. As such, we do not ex- 
pect the EROS survey to have detected the microlensing 
signal from our predicted sources; i.e. a faint population of 
SMC debris behind the LMC should not be of similar mag- 
nitude as the brightest LMC stars. It should be noted that 
the EROS team was the first to poin t out that most lensin g 
events were only seen for faint stars l|Tisserand et alj|2007l ). 

The OGLE II survey spanned a 4.7 deg 2 area centered 
on the bar region of the LMC. This is a somewhat smaller 
region than that covered by the MACHO survey. Following 
the MACHO strategy, blended sources were also considered. 
Over the 4 year span of the OGL E II survey (1996-2000), 2 
events were reported (LM C 1, 2, IWvrzvkowski et aT1l2009l : 
ICalchi Novati et alj|2009l , listed in Table 0). The follow-up 
OGLE III survey covered a larger 40 deg 2 area and reported 
2 new events (LMC 3, 4; listed in Table [2}, corresponding 
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to an optical depth of r = (0.43 ± 0.33) x 10 -7 , and 2 less 
proba ble events over a total sur vey time of 8 years (2001- 
2009) (|Wvrzvkowski et al.ll2011al ). The total event frequency 
is thus 4/12 ~ 0.33 events/yr. This frequency is smaller than 
that detected by MACHO, but the spatial distribution of the 
events is in agreement with the MACHO team's findings: the 
events are not concentrated solely in the bar region. The 
range of event durations is also similar. 

The OGLE III survey results imply a lower event fre- 
quency than observed by the MACHO team. This discrep- 
ancy might be explained by differences in the sensitivity 
and detection efficiency between the two surveys. The de- 
tection efficiency of the MACHO survey peaks at 50% for 
event durations of 100 days, whereas that of OGLE III peaks 
at 10-35% for the s ame event duration, depending on the 
density of the field (|Calchi Novati fe Mancinill201ll ). More 
significantly, the detection effici ency of the OGLE survey 
is lower in more crowded fields (|Wv rzvkowski et al.ll2011al ; 
ICalchi Novati fc MancinillioTlh : in contrast, the inner dens- 
est regions of the LMC disk is where the MACHO survey 
found the majority of their detections. Most of the LMC 
source stars are blue main sequence stars with little ex- 
tinction, as a result the MACHO blue and red bands have 
a much higher photon detection rate than the OGLE I- 
band. Another factor is that the MACHO team spent a 
larger fraction of its observing time on the MCs, whereas 
the OGLE team has higher sensitivity towards the Galactic 
bulge. Furthermore, the OGLE II and III survey threshold 
source magnitude was not as faint as that of the MACHO 
team (|Calchi Novati fc Mantinll201ll '). Note that, since the 
OGLE survey is more sensitive that the EROS survey, we 
opt to only directly compare the OGLE and MACHO results 
in this analysis. 

Note that in Tables [T]and we have not listed the sta- 
tistically corrected values of the blendecQ durations (Table 4, 
lAlcock et aLll2000l ). since the distribution of these durations 



is artifici ally narrowed, as pointed out by IBennettl (2005) 
(see also. lGreen fc Jedamzikll2002l ; lRahvarl2004h . This is fine 
for the calculation of the optical depth, but not for a study 
of the true distribution of durations. Instead, the durations 
listed in Table [T] are taken from Table 5 of lAlcock et alj 
(|2000h and Table 1 of iBennett et all (120051 ). The observed 
event durations reported by both teams range from (17-70) 
days, with similar average durations of 43 days for MA- 
CHO and 37.2 days for OGLE. No short duration events 
(lhr - 10 days) have been observed to date llAubourg et alj 
1 19951 ; iRenault et ail 1 19971 ; lAlcock et all 1 19961 . |2000D ruling 
out light objects (< lO^ 1 M©) as the prime constituents of 
the MW halo. 

Assuming the source stars are LMC disk stars, IBennettl 

(2005) recomputed the microlensing optical depth based on 
a refined list of candidates (Events 5 and 7 were also included 
in this reanalysis) and found r ~ 1.0 ± 0.3 x 10~ 7 . 

We aim here to explain the listed event duration, dis- 
tribution and event frequencies determined by the MACHO 
and OGLE teams (Tables [1] and [2j in a model where tidally 
stripped SMC stars exist behind the LMC's disk along our 



2 The term blending indicates that more than one star is 
enclosed within the seeing d isk of a given resolved object 
llCalchi Novati fc Mancinill201ll) 



line of sight. Note that a population of debris stars located 
further away from the LMC disk along our line of sight are 
not likely to have been been included in the EROS source 
selection of only the brightest (unblended) stars. The EROS 
null result is thus not a limiting factor to our analysis and 
we focus on the events detected by the MACHO and OGLE 
collaborations. 



3 THE SIMULATIONS 

We follow the methodology outlined in B12 to simulate the 
interaction history of the LMC-SMC-MW system using the 
smoo thed-particle hydrodynamics code, Gadget-3 (|Springell 
2005). We refer the reader to the B12 paper for details of the 
numerical simulations, but outline some of the more salient 
features here. 



3.1 Model Parameters 

Both the LMC and SMC are initially modeled with expo- 
nential gaseous and stellar disks and massive dark matter 
halos. The gaseous disk of the SMC is initially 3 times more 
extended than its stellar disk. The LMC has a total mass of 
1.8 x 10 11 M Q and the SMC has a total mass of 2 x 10 10 M 
(mass ratio of ~1:10). 

The interaction history of the LMC and SMC, indepen- 
dent of the MW, is followed in isolation for ~ 5 Gyr. The in- 
teracting pair is then captured by the MW and travels to its 
current Galactocentric location on an orbit defined by th e 
HST proper motion analysis of iKallivavalil et al.l (|2006al ). 
The MW is modeled as a static NFW potential with a total 
mass of 1.5 x 10 12 M©. In such a potential, backward orbital 
integration schemes show that the MCs are necessarily on 
their first infall towards our system, having entered within 
R200 = 220 kpc of the MW 1 Gyr ago. 

As mentioned in § 12.11 in this study we adopt the re- 
vised values of Vlsr fr om McMillan] (|201 lh and of V pec from 
ISchonrich et al' 

I (|201Ch. rather than the IAU standard val- 
ues as adopted in Kallivavalil et all (|2006al ) . The B12 models 
were designed to reproduce the 3D velocity vector and po - 
sition of the LMC as defined by IKallivavalil et all (|2006al ). 
thus we need to apply a velocity correction to the simulations 
in this study. In addition to the revisions to the Solar motion, 
our team has obtained a 3rd epoch of proper motion data 
using WFC3 on HST (Kallivayalil et al. 2012, in prep). With 
the new data, and better models for the internal kinematics 
of the LMC, the prop er motions have changed from those 
originally reported by IKallivavalil et all (|2006al lbh. Specif- 
ically, the total center of mass motion of the LMC is 56 
km/s lower than that reported by IKallivavalil et al. (2006a, 
378 km/s). 

To account for t hese updated valu e s, we s ystematically 
subtract out the old IKallivavalil et ail (|2006aT ) 3D velocity 
vector for the center of mass motion of the LMC from each 
particle in the B12 simulations and then add back the new 
Galactocentric 3D vector from Kallivayalil et al. (2012, in 
prep). 

If the LMC's total mass is 1.8 x 10 11 M©, as expected 



3 expe 

within the ACDM paradigm l|Bovlan-Kolchin et all [201lh . 
even this new lower velocity in a 1.5 x 10 12 M© MW model 
still represents a first infall scenario (Kallivayalil et al. 2012, 
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Table 1. MACHO Microlensing Events 



Event Event Duration 51 V~M h 
(days) 



1 17.3; 17.4 19.78 

4 39.5 21.33 

6 46.0 20.0 

8 33.2 20.20 

13 66.0 21.76 

14 c 53.3 19.48 

15 22 21.18 

18 37.9 19.55 

21 70.8 19.37 

25 42.7 19.04 



"Event durations arc taken from Table 1 of lBennett et alj i2005r i for events 4, 13 and 15, while the rest were taken from Table 5 (fits 

with blending) 

of lAlcock et al. I d200Ct) . rather than the statistically corrected 
durations used for the optical depth calculation. The average event duration 
is t e ~ 43 days and the event frequency is T ^1.75 events/yr over the 5.7 year survey. 
b V-band source star brightness from lBennett et alj d2005h ; lAkock et al.l <200ll, l2000l) 
c Event 14 is likely to be an LMC self-lensing event iAlcock et al.ll200lh 



Table 2. OGLE Microlensing Events 



Event 


Event Duration a 
(days) 


v M h 


Ol 


57.2 


20.65 ±0.013 


02 


23.8 


20.68 ± 0.014 


OS 


35.0 


19.52 ±0.1 


04 


32.8 


18.45 ±0.1 



"Event durations from OGLE (five p arameter fit): OGLE II resu lts, Ol and 02 from I Wvrzvkowski et aL (2009) 
and OGLE III results, 03 and 04 from I Wvrzvkowski etafl d2011al") . The average event duration is t e ~ 37.2 days and 

the event frequency is V ~ 0.33 events/yr. 
b Baseline V-band source star brightness from Wvrzvko wski et all |2009) for events Ol and 02 and I Wv rzvkowski et al. | [|2011al) for 03 

and 04 



in prep). Thus, the resulting dynamical picture has not 
changed by making this velocity shift. One can think of this 
velocity correction as changing the observer location, rather 
than the physics or dynamical picture. Moreover, similar re- 
visions must be applied t o the SMC's 3D velocity vector of 
iKallivavalil et ail j2006b) ; consequently we find that the rel- 
ative motion between the MCs is unchanged, at 128 ± 32 
km/s. In this study we are largely concerned with comput- 
ing quantities that are dependent on the relative motion 
of LMC stars (lenses) with respect to the background SMC 
stellar debris (sources). The exact value of the 3D velocity 
vector of the LMC is thus not a major source of error in this 
study. 

We caution that the B12 models do not reproduce ev- 
ery detail of the Magellanic System. In particular, the ex- 
act position and velocit y of the SMC as determined by 
IKallivavalil et all (|2006bl ) is not reproduced. The SMC's po- 
sition, however, does not affect the analysis in this study 
since we are concerned with the field of view centered on 



the LMC disk. Moreover, the goal of this study, and that of 
B12, is not to make detailed models of the Magellanic Sys- 
tem, but rather to illustrate a physical mechanism (namely 
tidal interactions between the MCs) that can simultaneously 
explain the observed large scale morphology, internal struc- 
ture and kinematics of the Magellanic System and, perhaps 
surprisingly, also the nature of the observed MACHO mi- 
crolensing events towards the LMC. 



3.2 Definition of Model 1 and Model 2 

In this study we examine the final state of the simulated 
Magellanic System for the two favored models for the orbital 
interaction history of the MCs from B12 (Model 1 and Model 
2). In B12 we characterized the gaseous debris associated 
with the tidal LMC-SMC interactions. Here we characterize 
the resulting stellar counterpart to this gaseous system. 

Model 1 and Model 2 differ in terms of the number 
of pericentric passages the SMC has made about the LMC 
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since they began interacting with each other ~6 Gyr ago. 
In Model 1, the SMC completes 3 passages about the LMC 
over the past 6 Gyr, but it never gets closer than ~20 kpc 
to the LMC. In Model 2 the SMC completes an additional 
passage about the LMC. This additional passage results in 
a direct collision between the MCs and the formation of a 
pronounced gaseous bridge connecting the MCs. A Bridge is 
reproduced in Model 1, but is of much lower density, having 
formed in the same tidal event that formed the Magellanic 
Stream. 

In B12, we showed that both of these models are 
able to reproduce the global large scale gaseous structure 
of the Magellanic System (Bridge, Stream and Leading 
Arm). Overall, however, Model 2 provides significantly bet- 
ter agreement with the internal structure and kinematics of 
the LMC. 

In particular, in Model 2, the LMC's bar (i.e. the pho- 
tometric center of the LMC) is offcenter from the geometric 
and kinematic center of the underlying disk and is warped 
by ~10 degrees relative to the disk plane, as observed. This 
is a consequence of a recent direct collision with the SMC. In 
contrast, in Model 1, the bar is centered and coplanar with 
the disk-plane. Model 2 is thus likely to provide a fairer as- 
sessment of the role of internal distribution of the LMC disk 
stars as lenses for the predicted microlensing event proper- 
ties. 



4 THE STELLAR COUNTERPART TO THE 
MAGELLANIC STREAM AND BRIDGE 



In lBesla et ail (|2010l ) and B12 we proposed that tidal forces 
from the LMC acting on the SMC are sufficient to strip 
out gaseous material that ultimately forms the Magellanic 
Stream without relying on repeated encounters with the 
MW. However, the dominant mechanism that removes mate- 
rial from the MCs is s till a subject of debate; stellar outflows 
l|Nidever et all 120081; [Ol ano 2004 ) and ram pressure strip- 
ping (|Mastropietro et all 120051 ; iMoore fc Davis! 1 19941 ) have 
also been proposed as possible explanations. Distinguishing 
between these formation scenarios is critical to the devel- 
opment of an accurate model for the orbital and interac- 
tion histories of the MCs and to our understanding of how 
baryons are removed from low mass systems. 

In the B12 model, LMC tides are also able to remove 
stars from the outskirts of the SMC disk in addition to the 
gas. In contrast, stars are not expected to be removed in 
the ram pressure or stellar outflow models. The detection 
of stellar debris in the Magellanic Stream would prove con- 
clusively that the Stream is in fact a tidal feature, ruling 
out models based on hydrodynamic processes. The theoret- 
ical prediction of a stellar counterpart to the Magellanic 
Stream, Bridge and Leading Arm is not novel; it is ex- 
pected in all tidal models f or the Magellanic Sy s tem (e.g., 
Gardiner fc Noguchj 1 1996; IConnors et al.l 120061 ; iLin et all 
19951 ; iRuzicka et ahlboiol ; biaz fc Bekkill2012l ). 

Recently, Weinberg (2012, in prep) found that by us- 
ing a kinetic theory approach to solve the evolution of gas, 
gas undergoing ram pressure stripping by some diffuse am- 
bient material is heated and ablated away, rather than being 
physically pushed out of the disk. As such, one would not 
expect ram pressure stripped material to be comprised pre- 



dominantly of neutral hydrogen. These preliminary results 
suggest that it is unlikely that ram pressure stripping by the 
ambient MW hot halo is responsible for the formation of the 
Magellanic Stream. This further implies that the location of 
the Stream on the plane of the sky cannot be used to pre- 
cisely constrain the past orbital trajectory of the MCs, since 
its orientation in a ti dal model will dev iate from the path 
followed by the MCs (|Besla et alj|2010l ). Note that, in the 
Appendix of B12 we also argue against a stellar outflow ori- 
gin based on the low metallicities observed along the Stream 
l|Fox et al.|[2010h . 

Ultimately, testing the tidal model has proven chal- 
lenging as previous searches for stars in sections of the 
Stream have yielded null resu lts |Guhathakurta fc Reitzell 
1 19981 ; iBrueck fc Hawkinj|l983l ) . We argue that this is likely 
because these studies were insufficiently sensiti v e. Bas ed on 
their null detection, iGuhathakurta fc Reitze] (|l998) con- 
strain the stellar-to-gas ratio to <0.1 in the high gas density 
clumps they observed in the region of the Stream known as 
MSIV (material between Galactic Latitudes of -70° and - 
80° a nd Galactic Longitudes of 90° and 45° ; iPutman et alj 
2003). Here we will compare the B12 model results with 
thes e constraints. 

Im 

ufioz et al] (|2006l ) claim to have detected LMC stars 
in fields near the Carina dSph galaxy, over 20 degrees away 
from the LMC. The identified stars have mean metallicities 
~1 dex higher than that of Carina and have colors and mag- 
nitudes consistent with the red clump of the LMC. They also 
argue that red giant star candidates in fields between the 
LMC and Carina exhibit a smooth velocity gradient, pro- 
viding further evidence of the kinematic association of these 
distant stars to the LMC. Such a velocity gradient would 
also be expected of tidal debris stripped from the SMC and 
the metallicities do not rule out an SMC origin (Ricardo 
Munoz, private comm unication, 2011). In line with this idea, 
iKunkel et al.l l|l997h argued that a polar ring of carbon 
stars from the SMC e xists in orbit about the LMC. Also, 
IWeinberg fc Nikolaevl (|200lh inferred the possible presence 
of tidal debris about the LMC based on Two-Micron All 
Sky (2MASS) data and suggest that this may have impor- 
tant implications for microlensing studies, depending on its 
spatial distribution. Such observations provide evidence for 
extended stellar distributions associated with the MCs that 
could potentially be interpreted as tidal debris. In this sec- 
tion we will attempt to reco ncile these observations with the 
tidal B12 model. Note that iMaiewski et al.l (|2009| s t have re- 
cently argued for the existence of an extended stellar halo 
around the LMC (David Nidever, private communication), 
but we do not model this explicitly. 

In B12 we noted that there is a significant transfer of 
stellar material from the SMC to the LMC expected in our 
proposed model. W e also showed that the recent detection by 
lOlsen et al.l (|201ll ) of a metal poor population of RGB stars 
in the LMC disk with distinct kinematics from the mean lo- 
cal velocity fields of the LMC is consistent with the proper- 
ties of the transferred stars in the B12 model. A marginally 
kinematically dist inct popula t ion of carbon stars has also 
been reported by iGraff et all (|200Ch , which they also sug- 
gest could be tidal debris in the foreground or background 
of the LMC that could account for the microlensing events 
observed toward the LMC. 

If stars are transferred from the SMC to the LMC there 
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must also be tidal debris in the bridge of gas that connect s 
the MCs, kno wn as the Magellanic Bridge (|Kerr et al.lll954h . 
iHarrisl l|2007h examined the stellar population in the Magel- 
lanic Bridge and was unable to locate an older stellar popu- 
lation, concluding that the detected young stars (< 300 Myr 
old) were formed in-situ in the gaseous bridge rather than 
being tidal debris s tripped from the SMC. More recently, 
iMonelli et alj (|201lh reported the detection of a population 
of stars located across the Bridge that could be compatible 
wi th the o l d dis k population of the LMC. The observations 
bv lHarris! (|2007T ) focused on the leading ridgeline (location 
of the highest gas density) of the Magellanic Bridge, which 
would experience the maximal ram pressure. It is possible 
that ram pressure from the motion of the MCs through 
the ambient halo gas of the MW may have displaced the 
gas in the bridge from the t idal stellar population. In fact, 
ICasetti-Dinescu et al.l (|2012f ) recently found a population 
of you n g star s offset by 1-2 degrees from the HI ridgeline. 
IHarrisl ^2007^ constrained the stellar density of a possible 
offset stellar population using 2MASS observations, but the 
2MASS sensitivity limit of 20 Ks mag/arcsec 2 is likely too 
low to have observed the expected stellar bridge. 

In the following we quantify the properties and observ- 
ability of the predicted stellar counterparts to the Stream, 
Bridge and Leading Arm Feature for the B12 model. 

4.1 The Extended Stellar Distribution of the 
Magellanic System 

Figures [1] and [2] show the Hammer- Aitoff projection of 
the stellar surface density and the corresponding Johnson 
V band surface brightness for Models 1 and 2, respectively. 
The expected Johnson V band emission is created from the 
simulations usin g the Monte Carlo dust radiative transfer 
code SUNRISE jjonsson fc Primaddl201(ih . The bulk of the 
stellar material originates from the SMC. To model this stel- 
lar debris we ad opt a Maraston stellar population model 
|Marastonll2005l ) and the oldest stars are assumed to have 
an age of 8.7 Gyr. SMC (LMC) stars are assigned metallic- 
ities of 0.2 (0.5) Zq, consistent with their observed values. 

The average V band signal in the Stream is > 30 
mag/arcsec 2 for both Models, well below existing observa- 
tional thresholds. This explains why the stellar counterpart 
has not been observed. The simulated stellar stream is more 
pronounced in Model 1 than in Model 2; this is because in 
Model 2, the Stream is older and the stars have more time 
to disperse. 

In particular, the stellar distribution in the Bridge 
region is well below 2MASS sensitivities (max of 20 
mag/a r csec 2 ) in both models, making it unsurprising that 
IHarrisl (|2007T l was unable to detect any offset stellar coun- 
terpart for the Bridge. The Bridge is clearly visible in the 
Vband in Model 2, but this is largely because of in-situ star 
formation (as the Bridge is known to be forming stars to- 
day). This is an important distinction between Models 1 and 
2; in Model 1 stars are not forming in the Bridge, whereas in 
Model 2 the Bridge is higher gas density feature as a result 
of the recent direct collision between the LMC and SMC. 
In Model 2, the Bridge is formed via a combination of tides 
and ram pressure stripping, owing to the passage of two 
gaseous disks through one-another. This explains why the 
stellar density of old stars in the Model 2 Bridge is low. No- 



tice that the stellar distribution in the Bridge is fairly broad: 
only the youngest stars exist in the high density "ridge" of 
the Bridge. The surface density of old stars in an y given 
sectio n of the bridge is low, making it plausible that IHarrisl 
|2007h would have been unable to identify any older stars in 
the small fields they sampled along the Bridge ridgeline. 

The two Models make very different predictions for the 
distribution of material leading the MCs. In Model 2, the 
leading stellar component is brighter than the Stream it- 
self. This leading component is not dynamically old, having 
formed over the past 1 Gyr and should therefore retain kine- 
matic similarities to the motion of the overall Magellanic 
System. 

To compa re against the upper limits on t he star-to-gas 
ratio placed by Guhatha kurta fc Reitzell (|l998T ) for the Mag- 
ellanic Stream, we consider a high gas density section of the 
Model 1 simulated stream. The chose n region would corre- 
spond to the MSI region as defined bv lPutman et all (|2003f) 
(between Magellanic Longitudes of -30 and -40). We choose 
to examine Model 1 because the gas and stellar densities in 
the modeled stream are higher than in Model 2: if Model 1 
is consistent with the limits then Model 2 should be as well. 

Figure [3] shows the simulated gas and stellar surface 
densities in the MSI region for Model 1 in Magellanic 
Stream coord i nates. This coordinate system is defined by 
iNidever et alj (|2008T ) such that the equator of this system 
is set by finding the great circle best fitting the Magellanic 
Stream. The pole of this great circle is at (l,b) — (188.5°, 
-7.5°) and the longitude scale is defined such that the LMC 
(l,b = 280.47°, -32.75°) is centered at Magellanic Longitude 
of 0°. From Figure [3] we find that the peak gas surface den- 
sity is ~0.3 Mq/pc 2 and the stellar surface density peaks 
at ~0.01Mq/pc 2 . The resulting stellar-to-gas ratio of 0.03 
is well within the observational limits. 

We thus conclude that the properties (density, surface 
brightness) of the predicted stellar counterpart of the Mag- 
ellanic Stream and Bridge from both Models 1 and 2 of the 
B12 study are not at odds with the null-detections from ex- 
isting surveys. 

4.2 The Source and Lens Populations 

In Figures 3] and [S] we illustrate the stellar distribution of 
the Magellanic System in Magellanic Coordinates, a vari- 
ation of th e galactic coordinat e system where the Stream 
is straight (|Nidever et all 120081 ). for both Models 1 and 2. 
The top panel shows the star particles from both of the LMC 
and SMC whereas the bottom panel depicts only stellar par- 
ticles from the SMC. The location of the LMC is indicated 
by the white dotted circle. There are clearly SMC stars in 
the same location as the LMC disk when projected on the 
plane of the sky. The distribution of stripped stars in Model 
2 differs markedly from that in Model 1; in Model 2 the 
SMC collides directly with the LMC, the smaller impact pa- 
rameter allows for the removal of stars from deeper within 
the SMC's potential. Consequently, the stellar debris in the 
Magellanic Bridge is more dispersed than in Model 1. 

In Figures|6]and[7]we zoom in on a 10 degree (~ 10 kpc) 
field centered on the LMC for Models 1 and 2, respectively. 
In the top panels we plot the stellar density of LMC particles 
in that field. These are the lenses for this analysis. Overplot- 
ted are the location of the MACHO and OGLE microlensing 



© 2012 RAS, MNRAS 000. rTH28l 



10 Besla et al. 





Figure 1. Hammer- Aitoff projection of the stellar density and Johnson V band surface brightness of the Magellanic System is plotted 
for Model 1. The orbital trajectory of the LMC(SMC) is indicated by the solid(dashed) white lines in the stellar density plots. The scale 
is the same in both images. 



candidates, labelled by their event number. In the middle 
panel we show the same field but plot only the SMC stellar 
particles; this is the surface density of sources. Again, the lo- 
cation of candidate microlensing events are marked and the 
white dotted ellipse indicates the extent of the LMC disk 
(this is not a circle in these images since the x and y axes 
are scaled differently than in the previous figures). There 
are significantly more sources in Model 2 than in Model 1: 
we correspondingly expect a higher microlensing event fre- 
quency in Model 2. 



Since the density of sources and lenses varies across the 
field, we have gridded the field and will compute quantities, 
such as the event frequency and duration, in each grid cell. 
The gridding and cell numbers are illustrated in the Figure [8j 
(plotted over the surface density of sources for Model 2) and 
will be referred to as marked throughout this text. Each grid 
cell is 2.5 x 3.75 degrees in size (at the distance of the LMC, 
1 degree ~ 1 kpc), representing the field of view over which 
we are computing the relevant microlensing quantities listed 
in § 12.11 The grid cells for comparison with the MACHO 
survey are indicated by the red box. The highlighted green 
region indicates cells used for comparison with the larger 
OGLE survey. 



5 THE PREDICTED MICROLENSING 

EVENTS IN THE B12 LMC-SMC TIDAL 
MODEL 

Next we compute the relevant quantities to determine the 
expected microlensing event durations and event frequencies 
for both Model 1 and 2 of B12, using the properties of the 
simulated source and lens populations. 



5.1 Effective Distance 

We begin our analysis by first computing the effective dis- 
tance between the identified lenses and sources in each grid 
cell, following equation |T|). The mean values of Ds, Dl and 
Dsl per grid cell are listed in Table [3] The weighted aver- 
age of D across the entire face of the LMC disk (weighted 
by the number of source stars in each grid cell) is D = 9.5 
(4.5) kpc for Model 1 (2). Sources are on average farther 
away from the lenses in Model 1 than in Model 2. Taking 
these average distances and Ei cns = IOOMq/pc 2 as the av- 
erage surface density in the central regions of the LMC, and 
following equation ((3), yields r ~ 7 x 10~ 7 for Model 1 
and t~3x 10~ 7 for Model 2. This suggests that the lens- 
ing probability is higher for Model 1, since on average the 
sources are located at a larger D. However, the observable 
of microlensing surveys is the number of events they detect 
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Model 2: Johnson Vbond 



0.8 2.0 




Figure 2. Same as Figure[T]but for Model 2. The stellar leading arm feature differs markedly from that of Model 1. The stellar density 
in the Bridge is also higher, however the majority of these stars are young, having formed in-situ in the high gas-density Bridge as a 
result of the recent SMC-LMC collision of Model 2. 




280 285 290 295 300 305 310 280 285 290 295 300 305 310 
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Figu re 3. The gase ous surface density (left) and stellar surface densit y (right) of a p ortion of the Magellanic Stream known as MSI 
llPutman et aljfeoOSh is plotted in Magellanic Stream Coordinates (Nidever et al. 2008t) . The simulated stellar stream is coincident with 
the gaseous stream and the stellar (gas) densities peak at a value of ~0.01 (0.3) Mq/pc 2 . 



per year. Obviously this number depends on the product of 
the number of sources and the optical depth. Since Model 
1 and Model 2 have different numbers of sources, the opti- 
cal depth alone does not tell the full story of which Model 
gives a larger number of microlensing events. Instead we will 
compute the event frequency explicitly. 



The SMC debris is not at a uniform distance from the 
LMC; consequently, there is a range of values for D per grid 
cell. The main grid cells where microlensing events have been 
observed are cells 7, 8, 9, 12, 13 and 14 (as highlighted by the 
red box in Figure In Figures 151 and IT01 we illustrate the 
probability distribution of D for SMC stellar debris (sources) 
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Figure 4. The stellar surface density of the Magellanic System is plotted for Model 1 in Magellanic Stream Coordinates. Magellanic 
Lo ngitude of 50 corresponds to a Galactic Longitude of Oo, and the LMC (l,b = 280.47°, -32.75°) is centered at Magellanic Longitude of 
0° llNidever et al . 2008, 2O10h. Lines of constant Galactic (l,b) are overplotted as dotted lines. The orbital trajectory of the LMC(SMC) 
is indicated by the solid(dashed) blue lines. The thick black line traces the true location of the Magellanic Stream on the plane of the 
sky. All particles are plotted in the top panel (this is the stellar counterpart of the top panel of Figure 7 in B12), whereas only particles 
originally part of the SMC are plotted in the bottom panel. There are clearly SMC stellar particles in the same field as the LMC disk 
(indicated by the thick black circle). 



relative to the mean distance of LMC lenses in each of these 
main cells. 



5.2 Relative Transverse Velocity of the Sources 

The event duration is dependent on the relative velocity of 
the source and lens perpendicular to our line-of-sight and 
projected in the plane of the lens (Vj_). We compute the 
relative velocity for all sources using equation (JHJ. 

There is a large spread in the possible relative velocities 
within each grid cell: the value for the lenses (Vl) is averaged 
per grid cell, but we keep track of the Vs for each source par- 
ticle individually in order to determine quantify this spread. 
In Figures [Til and [T2l we plot the normalized probability dis- 
tribution of V± for Models 1 and 2 in the central 6 grid 
cells where microlensing events have been detected. A wide 
range of relative velocities are possible, although velocities 
in excess of 300 km/s are unlikely. 

We integrate over the probability distribution of V± for 
each grid cell to determine the average (Vj_) per grid cell. 
These values are listed in Tables and 0] for Models 1 and 
2, respectively. The average weighted velocity across all grid 
cells for Model 1(2) is (V±) ~ 120(150) km/s; the kinematics 
are similar in both models. Given our massive LMC model 
(Mlmc = 1.8 x 10 11 Mq) and their relative spatial proximity 
to the LMC's center of mass (D < 10 kpc), such relative 



velocities imply that these debris stars are currently bound 
to the LMCq 

Fairly high relative velocities (60 — 200 km/s) are ex- 
pected in both models, but on average the relative velocities 
in Model 2 are higher than in Model 1 (150 km/s vs 120 
km/s). In Model 1, source stars represent tidally stripped 
stars that are in orbit about the LMC, forming a coher- 
ent, thin arc of stars behind the LMC. In Model 2, most 
of the source stars were removed during a violent collision 
between the two galaxies. This low impact parameter, high 
speed encounter allowed the LMC to strip stars from deep 
within the SMC's potential, naturally explaining the higher 
relative velocities of the stellar debris in Model 2. 



5.3 Event Duration 

Following equation ([9]), we compute here the average event 
duration ({t e }) per grid cell. As discussed in the previous sec- 
tions, the stripped SMC stellar particles (sources) exhibit a 
wide range of positions and velocities. Consequently, they 
will also exhibit a range of event durations. We thus com- 
pute a normalized probability distribution for i e , based on 

3 The escape speed at 20 (10) kpc from the LMC in the B12 
model is 220 (276) km/s 
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Figure 5. Same as Figure [4] except for Model 2. The top panel is the stellar counterpart of the bottom panel of Figure 7 in B12. 
During the collision between the LMC and SMC stars are stripped from the SMC from deeper within its potential than in Model 1. The 
corresponding debris has higher energy and appears more dispersed spatially than in Model 1. 



the probability distribution of the quantity D ' 5 /V±_ for the 
stripped SMC stellar particles in each grid cell. We plot the 
distributions for the central 6 cells of interest in Figure [T3] 
for Model 1, and Figure [Til for Model 2. 

Typical event durations range between 20-150 days for 
Model 1 and 20-100 days for Model 2. Short duration events 
(< 15 days) are not expected in these models, consistent 
with observations. Both models predict event durations in 
excess of 70 days, the longest event duration observed; how- 
ever event durations in excess of 100-150 days are rare. 

The observed spread in durations has often been at- 
tributed to the lensing mass function. Here we have assumed 
a single average value of (M^ s } for the lens and still obtain 
a large variation in the event duration in a given field of 
view. This is because of the large spread in relative speeds 
and effective distances of the tidal debris that represents the 
source population. The spread in the modeled durations per 
grid cell also explains why events such as MACHO Event 
13 and Event 1, which represent the longest and shortest 
events, respectively, can be closely located spatially. 

Integrating over these probability distributions yields 
an estimate for the average event duration, (i e ) (equation [9}, 
per grid cell, as listed in Tables 0and |5]for Models 1 and 2, 
respectively. Overall the event durations predicted by Model 
1 are longer than those of Model 2. 

The weighted average of the simulated event durations 
over the entire face of the disk yields: 80 days for Model 1 and 
42 days for Model 2. The shorter average event durations 
of Model 2 are more consistent with both the OGLE and 



MACHO results and are a direct consequence of the higher 
relative speeds determined in Model 2 (§ I5,2[I , 



5.4 Event Frequency 

Following equation fl6)l . we compute the average event fre- 
quency, (T) , per grid cell. We first compute normalized prob- 
ability distributions for F by accounting for the variation in 
the quantity D°' 5 V± exhibited by the stripped SMC stellar 
particles in each grid cell. These are plotted for the central 
6 grid cells in Figures [Tol and [TBI 

Integrating over these probability distributions yields 
the average event frequency in each grid cell; these are listed 
in Tables [3] and [4] for Models 1 and 2, respectively. 

Assuming 100% detection efficiency, the total predicted 
event frequency in Model 1 is T ~ 0.37 yr -1 for the entire 
LMC disk and 0.24 yr -1 for the central 6 grid cells repre- 
senting the area covered by the MACHO survey (red region 
in Figure [HJ. These values are much lower than expected for 
the MACHO events (1.75 yr" 1 ). To compare to the OGLE 
survey, we examine the cells in the larger green region high- 
lighted in Figure[S] We find an event frequency of 0.29 yr -1 , 
consistent with the results from the OGLE survey, which 
found 4 events over 12 years, implying an event frequency 
of 0.33 yr" 1 . 

For Model 2, the total event frequency across the entire 
face of the disk is ~5.17 yr -1 , 3.10 yr -1 for the central 6 
cells covered by the MACHO survey and 3.83 yr -1 for the 
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Figure 6. The stellar surface density of lenses (LMC stars) for Model 1 is plotted in the top panel for a 10 degree region centered on 
the LMC. The middle panel shows the surface density of sources (stripped SMC stars) located in the same field. The location of the 
MACHO and OGLE candidate microlensing events are marked by event number. 



OGLE survey region. These values are consistent with the 
MACHO survey result, but much higher than that of the 
OGLE survey. 



5.5 Detection Efficiency and Discrepancies 

between the OGLE and MACHO surveys 

So far we have not accounted explicitly for the detection 
efficiency of the MACHO and OGLE surveys in our deter- 
mination of the expected event frequencies. 

For the average simulated durations of ~ 20- 150 days for 
Models 1 and 2, the detection efficiency of the MACHO sur- 
vey is roughly 50%. As such, the observed event frequency 
should be roughly half that determined in the previous sec- 
tion for the central 6 grid cells (which well represent the 
MACHO survey area). That is, V ~1.54 yr" 1 for Model 2, 
which is in fact roughly the observed MACHO event fre- 
quency. Even taking a lower efficiency of 30% yields a con- 
sistent result of r ~1 yr" 1 . Since some events may turn out 
to be LMC self-len sing events (such as, MACHO event 14; 
lAlcock et al.|[200lh . this lower event frequency may be more 
accurate. However, the Model 1 event frequency is still an 
order of magnitude too low to match the MACHO observa- 
tions (0.07-0.12 yr" 1 ). 

For the OGLE survey, the maximal detection efficiency 



for the simulated average range of durations is 10-35%. The 
survey area for OGLE is also larger, so in the previous sec- 
tion we considered the event frequency expected across the 
region highlighted in green in Figure [8] finding a value of 
0.37 yr" 1 for Model 1 and 3.83 yr _1 for Model 2. Account- 
ing for the OGLE survey detection efficiency yields a sim- 
ulated event frequency of (0.03-0.10) yr -1 for Model 1 and 
(0.38-1.34) yr" 1 for Model 2. 

Generally speaking, microlensing should occur with the 
same probability on every type of source star. However, se- 
lection criteria adopted by the MACHO and OGLE imply a 
bias against detecting microlensing events on faint stars. In 
particular, the OGLE survey is less sensitive than the MA- 
CHO survey. Since our sources are expected to be faint SMC 
debris, a detection efficiency of 10% for the OGLE survey 
is likely appropriate. The observed event frequency by the 
OGLE survey of 0.33 yr - is consistent with the Model 2 
results of 0.38 yr -1 at such an efficiency. 

Discrepancies between the OGLE and MACHO surveys 
thus appear to be attributable to differences in the detection 
efficiencies of these surveys: specifically, their sensitivities to 
faint sources. As discussed in § 12.21 the EROS team non- 
detections towards the LMC is likely a result of their ap- 
proach of limiting their sample to only bright stars to avoid 
uncertainties owing to blended sources. 
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Figure 7. Same as Figure [6] except for Model 2. There are significantly more sources behind the LMC in Model 2 than in Model f. 
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Figure 8. The surface density of sources for Model 2 is illustrated with a grid ovcrplotted. Quantities will be computed within grid cells; 
each grid cell is assigned a number as indicated and will be referred to throughout the text. The red box indicates the cells (field of view) 
relevant for comparison with the MACHO survey. The green box highlights cells relevant for comparison with the larger OGLE survey. 



© 2012 RAS, MNRAS OOO.mf28l 



16 Besla et al. 



Model 



Cell 
Cell 



2 4 6 8 10 12 14 
D [kpc] 



Model 



Cell 12 



2 4 6 8 10 12 14 
D [kpc] 



Model 



Cell 13 
Cell 14 



m . n-i H- 1 . i- 




2 4 6 8 10 12 14 
D [kpc] 



Figure 9. The probability distribution of effective distances between individual stellar particles stripped from the SMC (sources) and the 
average lens position (LMC stellar disk particles) for individual grid cells, as indicated, for Model 1. Histograms have been normalized 
by the total number of stripped SMC stellar particles in each grid cell (listed in Table [3j ; integrating over the histograms yields a value 
of 1. For Model 1, Cell 7 does not have enough sources to compute a distribution. 




Figure 10. Same as Figure [9] except for Model 2. 



These above results favor Model 2, wherein the debris 
of SMC stars are non-uniformly dispersed behind the LMC 
disk, having been pulled out of the SMC or formed in-situ 
in gas stripped from the SMC after a collision between these 
two galaxies 100 — 300 Myr ago. 
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Figure 11. Probability distribution of the relative transverse velocity (Vj_) between individual stellar particles stripped from the SMC 
(sources) and the average velocity of LMC disk stars (lenses) for individual grid cells in Model 1. Histograms are normalized as described 
in Figure [9] 
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Figure 12. Same as in Figure llll but for Model 2. 
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Figure 13. Normalized probability distribution of the event durations (t e ) for selected grid cells in Model 1. 
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Figure 14. Same as Figure H"3l except for Model 2. 
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Figure 15. Normalized probability distributions of the event frequencies for selected grid cells in Model 1. 
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Figure 16. Same as Figure \TEl except for Model 2. Note, however, that the x axis is significantly larger than for Model 1. 
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Table 3. Microlensing Properties: Model 1 



Grid Cell 




( d l) 


( D ) 




<te) 


iVsource 


( r > 




(kpc) 


(kpc) 


(kpc) 


(Km / s ) 


(days) 


(10 4 ) 


(y r ~ ) 


it 




45.1 












at 




46.2 












3 




47.4 












4 


56.6 


48.7 


6.9 


435 


18 


0.25 


0.001 


5 


53.6 


49.8 


3.7 


115 


80 


6.0 


0.001 


et 




44.0 











0.0 


7*t 


61.0 


44.9 


11.9 


136 


76 


0.13 


0.001 


8*t 


57.9 


45.9 


9.4 


66 


108 


2.52 


0.020 


9 *t 


56.2 


47.1 


7.6 


50 


146 


5.66 


0.012 


1(1 


55.3 


48.7 


5.9 


64 


98 


8.68 


0.005 


lit 


59.4 


43.1 


11.7 


112 


40 


2.26 


0.007 


ia*t 


56.0 


44.0 


9.1 


147 


51 


4.53 


0.036 


13*t 


58.6 


45.1 


10.3 


132 


65 


7.81 


0.133 


14* t 


57.0 


46.2 


8.4 


83 


103 


9.06 


0.041 


15+ 


57.2 


47.6 


7.9 


(il 


138 


8.06 


0.010 


16 


59.9 


42.2 


12.4 


204 


47 


7.05 


0.008 


17 


59. 7 


43.2 


11.9 


183 


51 


11.70 


0.035 


18 


59.6 


44.1 


11.4 


169 


5a 


5.79 


0.030 


19+ 


58.7 


45.2 


10.3 


134 


53 


7.05 


0.025 


20 + 


57.7 


46.8 


8.8 


95 


80 


5.66 


0.007 


WEIGHTED AVERAGE 


57.9 


45.5 


9.5 


120 


79 


total 92.24 


total 0.37 


* MACHO Cells 


51.1 


40.9 


8.0 


88 


84 


total 33.22 


total 0.24 


t OGLE Cells 


57.6 


45.9 


9.1 


98 


90 


total 52.73 


total 0.29 


Note. — Asterisks indicate the 6 grid cells that 


cover reg 


ions where 


microlensing 


; events have 


been detected. 



Values for these six central cells (relevant area for the MACHO survey) are listed in the row marked '* MACHO 
Cells'. ' mark grid cells that cover regions relevant to the OGLE survey. The final row marked tf OGLE Cells' 
computes the total/average values relevant for the OGLE survey. Averages arc weighted by the number of source 
stars in each cell. Dashes indicate cells with no source stars. Event frequencies (T) quoted here assume detection 
efficiencies of 100%. 



Table 4. Microlensing Properties: Model 2 



Grid Cell 


( D s) 


( d l) 


( D > 


(VI) 


(t e ) 


N B0 urce 


<r> 




(kpc) 


(kpc) 


(kpc) 


{ km / s J 


(days) 


(10 4 ) 


(y r_ ) 


it 


57.8 


46.3 


9.2 


38 


28 


0.50 


0.001 


2+ 


50.5 


47.4 


2.9 


147 


23 


3.15 


0.010 


3 


51.0 


48.4 


2. (i 


130 


47 


7.81 


0.018 


4 


51.2 


49.3 


2.2 


69 


66 


4.28 


0.003 


5 


52.7 


50.1 


2.0 


144 


33 


2.89 


0.001 


6+ 


51.9 


45.3 


5.7 
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50 


22.01 


0.061 


7*t 


52.1 


46.2 


5.2 
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43 


23.40 


0.136 


g*t 


48.6 


47.1 


1.7 


66 


77 


88.82 


0.366 


g*t 


52.3 


48.5 


3.6 


88 


66 


13.46 


0.032 


10 


54.2 


49.2 


4,(i 


53 


96 


3.40 


0.003 


lit 


50.2 


44.7 


5.0 


177 


38 


758.10 


0.623 


12*t 


49.0 


45.4 


3.4 


157 


35 


669.59 


1.741 


13*t 


48.7 


46.6 


2.3 


123 


36 


159.40 


0.691 


14*t 


52.9 


47.6 


4.6 


130 


GO 


23.28 


0.133 


15t 


51.5 


48.3 


3.1 


129 


04 


7.43 


0.012 


l(i 


50.6 


44.8 


5.1 


158 


43 


26211.20 


0.776 


17 


49.8 


45.5 


4.0 


131 


47 


808.20 


0.372 


18 


50.4 


46.1 


3.9 


146 


41 


163.68 


0.157 


igt 


51.7 


46.9 


4.4 


181 


31 


28.32 


0.035 


20t 


49.7 


47.5 


2.3 


157 


24 


7.68 


0.002 


WEIGHTED AVERAGE 


50.2 


45.2 


4.5 


153 


42 


total 5416.52 


total 5.17 


* MACHO Cells 


49.1 


45.9 


3.1 


142 


40 


total 977.95 


total 3.10 


t OGLE Cells 


49.9 


45.6 


4.(1 


158 


39 


total 1797.72 


total 3.83 



Note. — Asterisks indicate the 6 grid cells that cover regions where microlensing events have been detected. 
Values for these six central cells are listed in the row marked as '* MACHO Cells'. The final row marked ( * OGLE 
Cells' computes the total/average values relevant for the OGLE survey. Averages are weighted by the number of 
source stars in each cell. Event frequencies (T) quoted here assume detection efficiencies of 100%. 
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6 OBSERVATIONAL CONSEQUENCES 

We do not claim in this work that all observed microlens- 
ing events have sources in an SMC debris population, but 
that such a population would be expected to have event 
frequencies and durations that are consistent with those 
observed. It is likely tha t the lenses do n ot all belong 
to the same population (|jetzer et al.l 120021 ') . In particu- 
lar, there is probably some con tribution to t he m icrolens- 
ing events from self-lensing: the iGvuk et al.l (|2000T ) disk + 
bar self-lensing models p r edict a contribution of 20% and 
ICalchi Novati fc Manciml <|201ll) were able to explain the 
OGLE events in a self-lensing model. Furthermore, the to- 
tal mass of the LMC's stellar halo is not well constrained; 
it may also contribute to the observed microlensing event 
frequency. In particular, on a first infall scenario the LMC 
is not expected to be truncated to the tidal radius, meaning 
that any putative stellar halo may be much more extended 
than originally considered. This is in line with the possible 
detection of LMC sta rs as far as 20 degr ees away from the 
LMC center of mass (|Munoz et al.l 120061 ). We further note 
that many of the initially reported microlensing events were 
subsequently determined to be binary lensing events; i.e. to 
have origins within the LMC disk. Thus, it is not unreason- 
able to consider a fraction of the lens populations as being 
within the LMC disk itself (|Sahull2003D . 

Here we outline observationally testable predictions 
from our model. 



6.1 Detecting the Stellar Counterpart to the 
Stream and Bridge 

The lBesla et al.l (|2010T ) and B12 models for the Magellanic 
System predicts the existence of a stellar counterpart to the 
Magellanic Stream and Bridge. We have demonstrated here 
that the surface brightness of the predicted stellar stream is 
too low to have been detected by surveys for Stream stars 
conducted to date (§HJ). 

There are a number of on-going optical searches for stars 
in the Stream, in cluding the NOAO O uter Limits Surve y 
jSaha et al.ll2010T l and the VMC survey (|Cioni et al.ll201lf ). 
that stand a good chance of detecting the predicted stel- 
lar stream. The NOAO Outer Limits Survey fields include 
non-contiguous pointings aligned with the gaseous Magel- 
lanic Stream and also sight lines in several directions from 
the MCs. For the published fields due north of the LMC, 
jSaha et all I2010T ) report the detection of LMC stars out 
past 10 disk scale lengths. The sensitivity of these data, 
corresponding to a surface brightness level in the Vband 
of -~32.5 mag/arcsec 2 (Abi Saha, private communication, 
2012), may be sufficiently sensitive to detect a spatially 
aligned stellar counterpart in their Stream fields. However, 
this survey will be unable to detect an offset stellar stream 
if one exists, as the search area is limited to regions of the 
highest gas density. The VMC survey has only one field on 
the Stream and one off the Stream. Future searches for stars 
would be more robust if designed to search contiguously 
across the Stream to detect changes in the magnitude of 
diffuse emission. 

The stellar Leading Arm in Model 2 (which is the model 
that best matches the kinematics and structure of the MCs) 
is predicted to be brighter than the trailing stream and is 



therefore likely the best location to look for tidal debris. Fur- 
thermore, a detection of stars in the Leading Arm may serve 
as a way to distinguish between Models 1 and 2. The stellar 
Leading Arm is expected to be offset from the gaseous Lead- 
ing Arm owing to the effect of ram pressure. Thus, searches 
for stars in the Leading Arm must survey a wide area of the 
sky surrounding the gaseous components. 

Overall, we find that the simulated stellar stream 
roughly follows the location of the simulated gaseous stream: 
no pronounced offsets exist in projection on the plane of 
the sky in simulations tha t do not include ram pressure/gas 
drag. The tidal models of iDiaz fc Bekkil |2012T l also suggest 
that a stellar counterpart must exist; however, they find that 
the tip of the predicted stellar stream is in fact offset from 
the observed gaseous Stre am, which could also exp l ain th e 
non-detection reported by iGuhathakurta fc Reitzell (|l99Sf ). 
Furthermore, if ram pressure effects were accounted for, the 
stellar stream would still appear as modeled here, while 
the gaseous stream may change location. As such, it ap- 
pears reasonable to search for the stellar counterpart of 
the Magellanic Stream by looking in regions displaced from 
the gaseous Stream by a few degrees, rather than just 
where the gas column den sity is highest. For example, 
ICasetti-Dinescu et al.l (|2012T ) find a population of stars that 
is offset by 1-2 degrees from the high-density HI ridge in the 
Bridge. 

From the toy model treatment of ram pressure stripping 
presented in the appendix of B12, we find that gas drag from 
the motion of the system through the MW's ambient gaseous 
halo does strongly modify the appearance of the gaseous 
Leading Arm. This will likely result in a severe discrepancy 
between the location of the gaseous and stellar Leading Arm 
features. It is therefore conceivable that stellar debris from 
the SMC may reach the Carina dSph fields ev en though the 
gaseou s Leading Arm does not, as suggested bv lMunoz et alj 
(l2006h . 

While ram pressure effects may also cause differences in 
the spatial distribution of stars in the Bridge relative to the 
gas, perhaps more significantly, in Model 2 the stellar debris 
in the simulated bridge is dispersed over an area ~ twice as 
large as that in Model 1. This is a consequence of the small 
impact parameter collision between the Clouds, which allows 
for the removal of stars from deep within the SMC potential 
well. In contrast, the gas is not similarly distributed, owing 
to hydrodynamic effects as the two gas disks pass through 
each other (e.g. gas drag). A spatially extended, low surface 
density stellar bridge is a generic prediction of an LMC-SMC 
collision scenario. 

6.2 Nature of the Source Population: SMC Debris 

We make a number of specific predictions for the proper- 
ties of the source stars of the observed microlensing events 
towards the LMC. First, these source stars originated from 
the SMC and should therefore be metal poor. This could be 
testable if spectra are obtained for the source stars. 

Second, the majority of these source stars should be 
old (age >300 Myr, roughly the age of the Bridge). Ages as 
high as 1 Gyr are not implausible based on isochrone com- 
parisons to the source population from the MACHO events 
(Tim Axelrod, private communication 2012). 

Finally, we also predict that the sources should be lo- 
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cated behind the LMC disk. A background source popula- 
tions has the advantage of b eing difficult to detect obser- 
vationally. IZhaol | 2000, f999) suggest that this population 
can b e identified by their expected reddening. IZhao et all 
|2000|)_ estimate this add itional reddening should be ~0.13 
mag. [Nelson et all l|2009h tested this expectation by compar- 
ing the color-magnitude diagrams (CMDs) of the observed 
sources to those of the LMC and detected no significant 
reddening. They consequently ruled out the possibility that 
these source stars may be behind the LMC disk. We note 
two caveats to this analysis. First, we consider here a source 
population originating from the SMC, which will have differ- 
ent CMDs than LMC stars. Second, new reddening maps of 
the LMC based on re d clump stars and R R Lyrae from the 
OGLE III data set bv lHaschke et all l|201ll ) reveal very little 
reddening in the bar region of the LMC. Since most events 
are concentrated in this region, we do not expect the source 
population t o exhib it significant reddening. Moreover, the 
iNelson et all (|2009h analysis only excludes a toy model in 
which the sources for all lensing events are in the background 
of the LMC and in which the extinction is uniform. For a 
more realistic model with a mixture of lens and source popu- 
lations and patchy extinction, the constraint is even weaker. 
We conclude that the existence of a source population of 
SMC debris located behind the LMC cannot be ruled out 
by the currently available data. 



6.3 Kinematics of the Source Population 

If the sources stars are indeed a population of tidally 
stripped stars from the SMC, we expect the majority of 
these stars to have distinct kinematics from the mean veloc- 
ities in the LMC disk. As discussed in B12, the line-of-sight 
velocities of these debris stars are, in some cases, discrepant 
from the local l ine-of-sight velociti es of the LMC disk. The 
observations of lOlsen et al.l (|201 lh provide suggestive evi- 
dence that such a metal poor debris population has already 
been observed. 

The determination of the proper motions of the sources 
may be a direct way to test our proposed model, as they 
should be discrepant from the mean proper motion of lo- 
cal LMC disk stars. In BI2 we have shown that the line- 
of-sig ht velocities of t his debris field are consistent with 



the lOlsen et al.l (|201lh results for both Models 1 and 2. 
lOlsen et alj l|201ll ) suggested two possible configurations for 
this debris that could explain the line-of-sight velocities: 1) 
a nearly coplanar counter-rotating disk population; and 2) 
a disk inclined from the plane of the sky at ~-19 degrees. 
Given that we expect the debris to exist behind the LMC 
disk, our models are more consistent with the second pro- 

pos ed configuration. 

ICasetti-Dinescu et al. (|2012l) recent ly measured the 
proper motions of the Idsen et all l|201ll ) data sample and 
concur that this population likely represents stripped SMC 
stars that are kinematically distinct from the LMC disk 
stars. Their analysis disfavors the coplanar counter-rotating 
disk scenario, but are consistent with an inclined orbital 
plane. They further find relative proper motions (i.e. tan- 
gential motions of the debris population relative to the LMC 
disk stars) of on average Afi a — —0.51 ± 0.30 mas/yr and 
Ajj,s = —0.03 ± 0.29 mas/yr. In this section we will test 



the predicted proper motions from our models against these 
results. 

Proper motions of source stars have been measured in 
a few cases, but most of thos e events are identifi ed as hav- 
ing MW lenses (e.g. Event 5; iDrake et alj|2004 ). An event 
of relevance is Event 9, a binary caustic crossing event 
jBennett et all 1 19961 ; iKerins fc Evans! 1 19991 ). The lens for 
this event is most likely a star near the LMC. This event was 
not included in our event frequency and dur ation analysis 
as it i s not included in the [Bennett] (|2005l ) and lAlcock et all 
(2000) conservative samples. The V± for the identified source 
(an A7-8 main sequence star) is low (~20 km/s). While this 
disagrees with the average V± velocities, we note that the 
velocities in cell 8 of Model 2 are lower on average; the dis- 
tribution in velocities in this cells peaks at V± ~ 66 km/s. 

We translated the distribution of relative transverse ve- 
locities (V±) between the sources and lenses in each grid cell 
(§ 15. 2p to a relative proper motion using the distribution of 
distances to the sources (§ 15. 1[) . We integrate over the result- 
ing probability distributions to determine the average proper 
motions per grid cell and plot the spatial distribution of 
these averages across the face of the LMC disk in Figure [T71 
Note that we plot only the magnitude of the proper mo- 
tion rather than the vector, as the vector is strongly model 
dependent. The weighted average of the simulated proper 
motions across all grid cells covering the face of the LMC 
disk is 0.44 mas/yr for Model 1 and 0.63 mas/yr for Model 
2. These values are consistent with th e average magnitude of 
the pr oper motion vector reported bv lCasetti-Dinescu et al] 
jjojg ) (~0.51 mas/yr). 

Follow up studies of the proper motions of the observed 
microlensed sources would be important tests to see if the 
kinematic properties of the sourc es are in fac t simil ar to 
the detected debris populat ion of lOlsen et al.l (|201ll ) and 
ICasetti-Dinescu et al.l |2012l ) and ultimately with our mod- 
els. The first epoch HST observations of all the MACHO 
sources have been taken in the 1990's, providing a signifi- 
cant time baseline. 



6.4 Spatial Distribution of Events and Event 
Durations 

From the results listed in Tables [3] and U it is clear that our 
models predict higher event frequencies in the central 6 grid 
cells that cover the LMC bar region. However, it is also clear 
that events are predicted to occur outside the bar region. 
The event frequency could be higher in these outer regions 
if the underlying source population is clumpy; i.e., such that 
there are high densities of source stars in fields not spatially 
coincident with the central bar (e.g. grid cell 16 in Model 2). 
Such o ff-bar events are not expe c ted in self-lensing models 
(e.g., ICalchi Novati et al.l 120091 ; ICalchi Novati fc Mancinil 
1201 ll ). which rely on the high surface density of lenses in 
the bar region and the warped nature of the bar itself. 

In § 15.31 we computed the average event durations in 
each of our denoted grid cells (listed in Tables [3] and |4| . In 
Figure [TH] we plot the spatial distribution of (t e ) across the 
face of the LMC disk. We find that both Models 1 and 2 
exhibit a trend such that, on average, event durations are 
longer at larger Magellanic Longitudes and Latitudes (i.e. 
North- West portion of the LMC disk). This is because t e 
is inversely proportional to the relative transverse velocity 



© 2012 RAS, MNRAS OOQ.[TH28l 



Microlensing Towards the LMC 23 





10 




8 




6 


~U 




.-i 




o 


4 


_1 




u 




E 


2 


o 




"53 




CTl 

o 





5 






_o 
z 




-4 




10 
8 


<D 


6 


"D 








s; 
o 


4 


_1 




(J 




'c 


2 


o 




QJ 




Qi 
D 





2 






-2 




-4 





4 2 0-2-4-6 

Magellanic Longitude 

Figure 17. The average relative proper motion between the sources and lenses is color coded in each grid cell across the face of the 
LMC disk for Model 1 (top) and Model 2 (bottom). The cells are defined in Figure [8] 



between the source and lens, which is on average smaller in 
that section of the disk (see the distributions of proper mo- 
tions in Figure [17]), whereas (D) is roughly constant across 
the disk. 

(VI) varies across the face of the LMC disk largely be- 
cause the LMC disk is rotating in a known clock-wise fashion 
and the disk is viewed roughly face-on. As such, the disk ro- 
tation is a dominant contribution to the tangential relative 
motion between the sources and lenses. Since the kinematics 
of the source population is largely determined by the rela- 
tive motion between the LMC and SMC (which is similar 
in both Models since it is dependent on the HST proper 
motions), both Models 1 and 2 exhibit similar behavior. 

If more events are observed at larger distances from the 
bar it may be possible to observe this average difference in 
durations. The currently observed events are too close to the 
central disk regions for this effect to be obviously discernible. 



6.5 Surface Density and Mass of the Source 
Population 

Model 1 fails to reproduce the observed event frequency, 
underestimating the rate by an order of magnitude. In the 
previous section we have illustrated that the kinematics in 
both models are roughly similar. Thus, the discrepancies be- 
tween Models 1 and 2 likely reflect the significantly different 



surface densities and mass of stripped SMC stellar debris 
behind the LMC disk in each model. Specifically, the num- 
ber of sources are too small and too smoothly distributed in 
Model 1. 

In Model 1 only a modest amount of stars are stripped 
from the SMC to the region behind the LMC disk; the stellar 
debris comprises only 0.2 per cent of the current LMC disk 
mass. In contrast, the mass of stellar debris in Model 2 is 
1.5 per cent of the LMC disk mass. 

The failure of Model 1 provides a lower limit of 7 x 
10 6 M for the mass of SMC debris behind the LMC disk. 
The distribution of this material could be clumpy, but the 
average surface density of this debris must be larger than 
that of Model 1, i.e. E dc bris > 0.03 Mq pc" 2 . 

Uncertainties in the star formation prescription adopted 
can change the amount of stellar debris predicted in the 
simulations. Regardless of such uncertainties, the limits dis- 
cussed here are still realistic constraints on the stellar den- 
sities required to reproduce the observed microlensing prop- 
erties. 

lOlsen et al.li|201ll ) find a kinematically distinct popula- 
tion of stars in the LMC disk field that comprise 5 per cent 
of their sample. Such studies can provide important tests of 
the limits placed here and better constraints on the density 
and distribution of the stellar debris from the SMC. 
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6.6 RR Lyrae in the Magellanic Stream, Bridge 
and Leading Arm 

The existence of a stellar counterpart to the extended 
gaseous components of the Magellanic System may be in- 
ferred via the detection of RR Lyrae stars through e.g. the 
ongoing OGLE surveys. The variability of these stars and 
their reliability as distance indicators make them readily 
identifiable tracers of a faint stellar population that can be 
associated with Magellanic System. 

Given the stellar mass in the Bridge, Stream and Lead- 
ing Arm from the simulations presented in this study, we 
can make a simple estimate for the expected number of RR 
Lyrae. 

Using the OGL E III Catalog for Variable Stars, 
ISoszvnski et"aH feoich find 2445 RR Lyrae in the SMC. 1933 
of these pulsate in the fundamental mode (RRab). If we as- 
sume that the RR Lyrae are uniformly distributed over the 
face of the SMC, we can estimate the number of RR Lyrae 
per unit stellar mass in the SMC. 

The total stellar mass of the SMC is 3 x 10 s L x = 

3 x 10 8 Mq (adopting a mass to light ratio of IMq/Lq). 
This means that there are 



l R 8 L // ae = 6.4 x 10- 6 RRLyrae/ Mq (18) 



In the B12 model the Stream, Bridge and Leading Arm 
all originate from the SMC. Thus, we can use this conversion 
factor to estimate the corresponding number of RR Lyrae in 
each of these components. The results are listed in Table [5] 

Number counts of RR Lyrae in the Stream are expected 
to be low in both models (< 40). These stars would also 
be scattered across the ~100 degree extent of the Stream, 
making the number density of RR Lyrae undetectably low. 

RR Lyrae in the Leading Arm are too low in Model 
1. However, the stellar counterpart to the Arm in Model 2 
is more massive (see § 14. 1[) . As such, the RR Lyrae counts 
predicted by Model 2 may be detectable: 160 RR Lyrae are 
expected over the 70 degree extent of the Arm. 

Detecting RR Lyrae looks most promising for the Mag- 
ellanic Bridge in Model 2, where we expect ~250 RR Lyrae. 
The quoted values are estimates that invoke a number of as- 
sumptions, the corresponding error bars are expected to be 
large. However, they do suggest that searches for RR Lyrae 
should focus on the Magellanic Bridge connecting the MCs. 

If RR Lyrae are detected, reliable distances can be ob- 
tained to map out the 3D structure of the Magellanic Bridge. 
Such observations can significantly constrain models for the 
recent interaction history of the Magellanic Clouds and val- 
idate our proposed theoretical model for the origin of the 
microlensing events. 
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Table 5. RR Lyrae Number Counts in the Magellanic System 



Component Stellar Mass Number of RR Lyrae a 

(1O 6 M ) 



SMC b 300 1933 

Model 1 Stream 6 38 

Model 2 Stream 4 26 

Model 1 Leading Arm 5 32 

Model 2 Leading Arm 25 160 

Model 1 Bridge 2 13 

Model 2 Bridge 40 250 



"Number of expected RR Lyrae pulsating in the fundamental mode (RRab). 
Values are computed from the stellar mass using the conversion factor of 
6.4 X 10 -6 RRLyrae/ Mq. 
b Stellar mass from Stanimir ovic et al.l J2004T) . RR Lyrae number counts from [Soszvnski et al. I J2010T) . 



6.7 Metal Poor Globular Clusters in the MW 
Halo 

The SMC is known to harbor an anomalously low number 
of old globular clusters. Although most SMC cl usters may 
have formed in two main events, 2 and 8 Gyr ago (Rich ct al 



2000), there are also clusters at intermediate ages (|Piatti 



20111 ). This is not true of the LMC, where a pronounced 
age gap exists; clusters in the LMC are either as old as 



the universe or 1-3 Gyr o ld (Da Costa 1991; Ivan den B crgh 
199ll ; iGirardi et aT]|l995l ; lOlszewski et al.lll996l : iPiatti et al 



2002). However, the oldest SMC cluster is at the young and 
metal-rich extreme of t he LMC globular cluster distribution 
(|Piatti fc Geislerll2012T ); there are no analogs of galactic halo 
lobular clusters in the SMcQ, but they do exist in the LMC 
Da Costal fl99ll ). 

We suggest that the fundamental difference in the clus- 
ter age distribution in the LMC and SMC might be a result 
of their interaction history. 

The recent dramatic collision between the two galaxies 
predicted by Model 2 could have resulted in the destruc- 
tion of globular clusters in the SMC, explaining its over- 
all paucityQ Furthermore, older clusters might have been 
stripped by the action of LMC tides over time. Some could 
be accreted by the LMC, but stellar material is not trans- 
ferred as efficiently between the two galaxies as gaseous ma- 
terial. The stripped clusters would instead most likely reside 
in the galactic halo, tracing the extended Magellanic Stream 
and Leading Arm on the plane of the sky. 

It is thus possible that old globular clusters in the galac- 
tic halo that exist at Galactocentric distances of ~50-120 
kpc and are spatially coincident with the Magellanic Stream 
or Leading Arm might have once resided in the SMC. If such 
clusters exist, their kinematics could confirm an association 
with the Magellanic System. 



4 The oldest cluster in the SMC is NGC 1 21, which is 2 Gy r 
y ounger than t he old est clusters in the LMC llShara et alJll998T) 

5 ICIatt et al.l d2010h find no evidence of dissolution of clusters 
younger than 1 Gyr, but have no constraints for older clusters. 



7 CONCLUSIONS 

We have demonstrated that a population of tidally stripped 
stars from the SMC located ~4 -10 kpc behind a lensing 
population of LMC disk stars can naturally explain the ob- 
served event durations, event frequencies and spatial distri- 
bution of microlensing events from the MACHO and OGLE 
surveys. These results favor the B12 Model 2 scenario for 
the interaction history of the MCs, wherein the MCs are 
on their first infall towards the MW and the SMC has re- 
cently collided with the LMC, leading to a large number of 
sources distributed non-uniformly behind the LMC disk. We 
note that the models presented in this study were not engi- 
neered to address the microlensing puzzle, but were rather 
constructed to explain the structure and kinematics of the 
Magellanic System more generally. Despite this fact, Model 
2 seems to provide a natural explanation for the observed 
properties of the reported microlensing events. 

The models presented in this study are consistent with 
the upcoming 3rd epoch HST proper motion data (Kallivay- 
alil et al. 2012, in prep). Note that, since the microlensing 
observables are dependent on the relative motions of the 
sources and lenses, the exact 3D velocity vector adopted for 
the system is not a major source of error. Both Models 1 and 
2 use cosmologically motivated models for the dark matter 
content of the SMC and LMC and can ex plain the origin of 
the Magellanic Stream |Besla et al.l [201Cl|). However, in B12 
we showed that only the direct collision scenario of Model 
2 reproduces the internal kinematics and structure of both 
the LMC and SMC. Support for such an off-center, moder- 
ate !_to_lfigh4mTUria^^ has recently been proposed 
bv lCasetti-Dinescu et all (|2012T ) based on the configuration 
of OB candidates in the LMC. In this study we find that 
Model 2 also provides a better match to the observed event 
durations and event frequencies reported by both the OGLE 
and MACHO microlensing surveys. Model 1, on the other 
hand, fails to reproduce the event frequency. 

The lower event frequencies reported by the OGLE sur- 
vey appears to be naturally accounted for by the survey's 
lower detection efficiencies and sensitivity to faint sources, 
compared to the MACHO survey. Assuming a detection ef- 
ficiency of 30-50% we find event frequencies of F ~ 1-1. 54 
events/yr in the central regions of the LMC disk in Model 
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2; this is consistent with the observed rate for the MACHO 
survey of T ~1.75 yr _1 . If instead we assume a much lower 
detection efficiency of 10%, we find event an event frequency 
of F ~0.38 yr _1 across a larger region comparable to that 
covered by the OGLE survey; this is consistent with the 
OGLE survey result of V ~0.33 yr~ . Model 1 underpredicts 
the MACHO observed event frequency by a factor of ~10, 
although this predicted rate may still be consistent with the 
OGLE results for a high detection efficiency (~0.29 yr" 1 
across the entire face of the disk, 100% efficiency). 

Models 1 and 2 both reproduce the observed range of 
event durations (t e ~ 17-70 days) without appealing to vari- 
ations in lensing mass. Instead, we find that the simulated 
sources exhibit a range of effective distances and relative 
transverse velocities within a given field of view owing to 
their tidal origin. On average, event durations are shorter 
for Model 2, reflecting the higher relative velocities of the 
stellar debris removed during a high speed, low impact pa- 
rameter collision between the Clouds. 

The simulated lensing events are expected to occur pre- 
dominantly, but not exclusively, in fields of view in the cen- 
tral regions of the LMC disk. This is in sharp contrast to 
self-lensing models. Self-lensing models may be able to ex- 
plain t he lower event frequency determined by the OGLE 
surve y (|Calchi Novati fc Mancinill201ll ; ICalchi Novati et all 
2009), but they fail to explain the existence of events in 
fields off the LMC's stellar bar (such as MACHO events 21, 
18); there is a strong decrease of the LMC self-lensing rate 
with increasing distance from the LMC stellar bar. The ex- 
act event frequency expected in off-bar fields in our proposed 
model strongly depends on the spatial distribution/surface 
density of the unde r lying source debris field. Studies like 
that of lOlsen et al.l (|201ll ) may be able to constrain this 
unknown. 

Model 2 can thus reproduce the observed microlensing 
quantities reported by both the MACHO and OGLE sur- 
veys. The EROS team chose to limit their sample such that 
they monitored only the brightest stars in the LMC; we do 
not expect that this approach could detect the lensing signal 
from a population of faint SMC debris, potentially explain- 
ing the EROS survey non-detection towards the LMC. We 
do not compare directly to the microlensing optical depth 
reported by these three surveys, as this is a derived quan- 
tity that depends sensitively on what is assumed about the 
source population. These surveys have assumed that the 
sources represent a population of LMC stars, whereas we as- 
sume here that the sources represent SMC stellar debris. As 
such, our derived optical depth cannot be compared directly 
to the values quoted by the MACHO, EROS and OGLE 
teams. Instead we have computed the observable quantities 
(event duration and event frequency) directly. 

The presented scenario makes a number of testable pre- 
dictions. In particular, invoking a stellar debris field as the 
source population implies that there must exist a stellar 
counterpart to the gaseous Magellanic Stre am (see also, 
iGardiner fc Noguchilfl99l ; iDiaz fc Bekkll2012l ) . Searches for 
stars in the Stream have thus far yielded null results. In § 2] 
we have shown that the predicted stellar stream has a sur- 
face brightness in the Vband that is > 34 mag/arcsec 2 . This 
is well beyond the sensitivities of these surveys and so our 
model does not violate known observational constraints. We 
further predict number counts of RR Lyrae in the simu- 



lated Magellanic Bridge, Stream and Leading Arm; number 
counts in the Bridge and Leading Arm in Model 2 may be 
detectable. The variability of RR Lyrae make them easily 
identifiable and their reliability as distance indicators pro- 
vides an accurate method for confirming association with 
the Magellanic System. The detection of RR Lyrae may thus 
provide the most direct method of confirming the tidal na- 
ture of the extended gaseous components of the Magellanic 
System and validating the theoretical model presented in 
this study. 

Model 1 and Model 2 also differ in the predicted spatial 
distribution of debris in the Bridge. In the Model 2 SMC- 
LMC collision scenario, LMC tides are able to remove stars 
from the inner regions of the SMC potential. This debris will 
consequently be more energetic and spread out over a larger 
area than that predicted by Model 1 or that traced by the 
gas. 

This model predicts that the sources should originate 
from the SMC and that they should make up more than 
~0.2 per cent of the LMC's disk mass and be distributed 
such that the surface density is larger than 0.03 Mq pc -2 . 
This lower limit is derived from the failure of Model 1 to re- 
produce the observed event frequencies. The Model 2 sources 
make up 1.5 per cent of the LMC's current disk mass and 
naturally explain the observed microlensing events. Such a 
populatio n of low-me t allicit y stars may have already been 
detected; lOlsen et al.l l|201ll ) have identified a population of 
low metallicity stars with distinct kinematics from the LMC 
disk stars in the same field of view. These stars comprise 
5 per cent of their sample and we have already shown in 
B12 that their line-of-sight kinematics are consistent with 
expectations of our Models 1 and 2. In this study we have 
made further predictions for the expected proper motions 
of these source stars with respect to local LMC disk stars, 
which could be observationally testable (§ [6T3J . 

In addition, owing to the rotation of the LMC disk, we 
find that there is a gradient in the relative tangential ve- 
locities of the source population across the face of the LMC 
disk. Consequently, we also expect to see, on average, a gra- 
dient in the observed event durations of the corresponding 
microlensing events (since t e oc V£ , § I6.4|l . It may be pos- 
sible to observe this gradient if a number of microlensing 
events are detected at locations a few degrees away from 
either side of the LMC's stellar bar. 

This study provides a natural explanation for the long- 
standing puzzle concerning the origin of microlensing events 
towards the LMC, which, if unexplained by normal stellar 
populations, would imply the existence of exotic MACHO 
objects populating the MW halo. The success of Model 2 
in reproducing the observed lensing event properties sup- 
ports a scenario in which the LMC and SMC have recently 
(<300 Myr ago) collided directly, causing the removal of 
stars from deep within the SMC potential and the formation 
of the Magellanic Bridge. Moreover, the agreement between 
our simulations and the MACHO and OGLE survey data 
also suggests that there exists a stellar counterpart to the 
Magellanic Stream and Bridge that has yet to be discovered. 
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